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ABSTRACT 

A simple model of gas accretion in young galaxy disks suggests that 
fast turbulent motions can be driven by accretion energy for a time t acc ~ 
2 (e 5 G M 2 / ^V 3 )°'° where e is the fraction of the accretion energy going into 
disk turbulence, M and V are the galaxy mass and rotation speed, and £ is the 
accretion rate. After t acc , accretion is replaced by disk instabilities as a source 
of turbulence driving, and shortly after that, energetic feedback by young stars 
should become important. The star formation rate equilibrates at the accretion 
rate after 1 to 2 t acc , depending on the star formation efficiency per dynamical 
time. The fast turbulence that is observed in high redshift starburst disks is not 
likely to be driven by accretion because the initial t acc phase is over by the time 
the starburst is present. However, the high turbulent speeds that must have been 
present earlier, when the observed massive clumps first formed, could have been 
driven by accretion energy. The combined observations of a high relative veloc- 
ity dispersion in the gas of z ~ 2 clumpy galaxies and a gas mass comparable to 
the stellar mass suggests that either the star formation efficiency is fairly high, 
perhaps 10 x higher than in local galaxies, or the observed turbulence is powered 
by young stars. 



Subject headings: galaxies: evolution — galaxies: high-redshift — galaxies: starburst 
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1. Introduction 

Deep surveys have detected Milky- Way-size galaxy disks at redshift z ~ 2 that con- 
tribute a large fraction of the stellar mass density and have star formation rates of order 100 
M yr" 1 (Daddi et al. 2004, 2007; Forster-Schreiber et al. 2006; Erb et al. 2006a,b; van 
Dokkum et al. 2006). They appear to be highly turbulent with gas velocity dispersions of 
40-80 km s" 1 (Forster-Schreiber et al. 2006, 2009; Genzel et al. 2006, 2008; Cresci et al. 
2009). They form stars in a small number of clumps with masses of around the turbulent 
Jeans mass, which is ~ 10 8 M (Elmegreen & Elmegreen 2005; Bournaud et al. 2007, 2008). 
They are also gas-rich, with gas-to-star mass fractions of order 30-60% (Tacconi et al. 2009). 
These features are in contrast to present-day disk galaxies that are relatively quiescent by 
comparison, having turbulent velocities of order 10 km s _1 (Dib et al. 2006), star-forming 
complexes typically smaller than 10 6 M Q , and star formation rates of order 1-3 M & yr _1 . 

The presence of gas does not necessarily imply star formation. Certain conditions have 
to be met before gas can become unstable to form self-gravitating clouds (e.g., Toomre 1964; 
Goldreich & Lynden-Bell 1965), and other conditions are required before these clouds can 
cool enough to condense into stellar objects (e.g., Elmegreen & Parravano 1994; Schaye 2004; 
Krumholz, McKee, & Tumlinson 2009; Prochaska & Wolfe 2009). These conditions have the 
effect of delaying the onset of star formation in a young galaxy, allowing the gas mass to 
build up without attrition until it reaches a tipping point, after which the conditions change 
and star formation is fast. The observed star-formation rate of 100 M Q yr~ x corresponds to 
nearly a whole galaxy disk of gas, several xl0 10 M Q , being converted into stars in only a few 
orbit times, which is several xlO 8 yrs. If the gas accretion rate is comparable to this star 
formation rate, then the initial build up of a disk is fast and star formation has to be very 
efficient, with a shorter gas consumption time than in local galaxies by a factor of ~ 5. If 
the accretion rate is less than the star formation rate, then there was probably a prior phase 
when most of the gas assembled before the star formation began. 

The properties of galaxy disks during the initial accretion and starburst phase are not 
understood. The large masses of star-forming clumps suggest that the gas was highly tur- 
bulent before the clumps formed. In this case, the primary source of turbulent energy would 
seem to be the accretion itself (e.g., Elmegreen & Elmegreen 2005; Forster-Schreiber et al. 
2006; Genzel et al. 2008; Dekel et al. 2009; Khochfar & Silk 2009; Klessen & Hennebelle 
2009). After the disk becomes unstable, further energy can come from disk self-gravity 
(Burkert et al. 2009; Bournaud et al. 2009). Here we investigate the duration of the initial 
turbulent phase when most of the energy comes from cosmological accretion, and we con- 
sider the cooling that leads to instabilities and star formation. We write equations for energy, 
mass, and star formation in dimensionless form, and then normalize these equations in a way 
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that is independent of the gas accretion rate and the efficiency of conversion from galaxy 
potential energy to disk turbulent energy. The equations are in Section 2, the solutions are 
in Section 3, and a comparison with observations is in Section 4. 



We simplify a galaxy during the short phase of initial accretion by considering a constant 
gas accretion rate dM^^/dt = £ into a disk with fixed radius R within a fixed potential V 2 . 
The energy accretion rate is approximately £ multiplied by the depth of the potential well, 
V 2 , assuming the thermal speed in the accretion flow is small compared with V. We let 
e be the efficiency of conversion of accretion energy into turbulent energy, E, in the disk. 
If e < 1, then some of the potential energy from accretion is lost to shocks and thermal 
radiation before it drives disk turbulence. The rate of input of turbulent energy is e^V 2 . If 
dark matter accretes along with the gas, then V 2 will increase with time. This is considered 
in Section HI 

An accreting disk dissipates turbulent energy at a rate approximately equal to the 
current energy divided by the perpendicular crossing time (e.g., Mac Low et al. 1998; 
Stone, Ostriker & Gammie 1998). This crossing time is t pcrp = H/a for disk scale height 
H = CT 2 /7rG£ disk = a 2 R 2 /GM disk . The turbulent speed is denoted by a = (2£/3M gas ) ' 5 
and the accreted disk mass is Mdi S k = ^fdisk^- Thus the dissipation rate is Ea/H = 
G(1.5M gSuS E)°- 5 MdiskR~ 2 ■ Putting accretion and dissipation together, the time-dependent 
equation for turbulent energy is 



The last term is the energy lost by locking up interstellar matter with velocity dispersion 
a into stars at a rate S. This counts for energy lost from the gas phase; i.e., the total gas 
energy decreases from this effect at a rate equal to the rate of conversion of gas into stars 
multiplied by the energy content of that gas. 

Star formation converts disk gas into stars. We take the disk gas mass to be M gas and 
the disk star mass to be M stars , so that M gas + M stars = M disk . No gas or stars are assumed to 
be ejected from the disk. The star formation rate is taken to be Sq times the gas mass M gas 
multiplied by the growth rate of a gravitational instability, which is u = 7rGS gas /cr in the 
linear instability theory. The coefficient Sq is an efficiency of star formation in a dynamical 
time; it is observed to be So ~ 0.01 on a variety of scales (Kennicutt 1998; Krumholtz & Tan 



2. An Accretion-Dissipation Equation 




(1) 
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2007) reflecting the small fraction of cloud mass that is dense enough to form stars and the 
rapid disruption of star-forming cores. The total star formation rate in the galaxy is then 

GM* , — GM 2 £ 
S = *_= = iW L6x_J5. (2) 

Note that accretion followed by disk shrinkage at constant M gas causes the gas surface den- 
sity to increase and the radius R in this equation to decrease. Shrinkage is expected if 
gravitational instabilities drive turbulence, because then turbulent energy is taken from 
gravitational binding energy, and the dissipation of this energy requires settling of the gas 
to a smaller radius. In the following discussions, the length scales will be normalized to the 
galaxy radius R. Then disk shrinkage should be viewed as a larger effective value of Sq, 
according to equation (J2J). 

Because of star formation, the gas mass in the disk varies with time as 

^= { -S, ,3) 

and the star mass varies as 

dM sta , IB g ^ 



dt 



We now normalize the important physical parameters, writing time in units of R/V, gas 
or disk mass in units of galaxy mass, M , and energy in units of M V 2 . Geometric terms 
of order unity are ignored, as are variations throughout the disk. There are two parameters 
in this model, £ and e. We can absorb them into our new variables to write the equations 
independent of them. Then we get general solutions to the equations that are independent 
of all parameters. The gas accretion rate £ is dimensional, so we introduce the dimensionless 
accretion rate, A = £R/(VM ), and define: 

[2RV 2 _ E , _ tA°- 5 V 

~ V 3 ~G~ A^e^MoV 2 ' * ~ l^R ' 

M disk 

. /l/f' Mggs _ » ri -^stars 



""disk /in.S,n.9K J/ ' gas a n 5 n ok 71 r 1 ""stars 



a' 



A°- 5 e°- 25 M ' gas A 05 e°- 25 M ' stars A°- 5 e°- 25 M 
( 2E' \ 5 S SR 



e 05 V \3M'J ' £ AVMq' 



Then equations ©-(jlj) can be written 

dE' 



(W 1 - K* (e'KJ °' 5 - S ' E '/K™ ( 5 ) 
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S' = S [M' o /E')°' 5 (6) 



dM' 



gas 
disk 



dt' 

dM' 



1 (7) 



gas 



1-S' 



dt 

We also define the useful quantities 



dt' 

S>. (9) 



dM' stais 



Q' = QA°- 5 /e - 25 = 3 - 5 a'/M' gas 
Mj eans = M Jeans A°- 5 /(e L75 M ) = 1.5n(a') A /M' gas 
J = ue°- 25 R/{A°- 5 V) = {2/3)°- 5 {M' g Ja') 



Here Q = ko /(tvGTH) is the Toomre instability parameter for epicyclic frequency k = 2°' 5 V/R 
in a flat rotation curve and mass column density £ = M gas / (ttR 2 ). The unstable Jeans mass 
is Mj eans = cx 4 /(GS). The growth rate u> was given above. All primed quantities are 
independent of £ and e. Note also that oj't' = ut. 

Equations ©-© were solved numerically. The initial conditions are M' gas = 0, M' disk = 
0, and E' = 0. The solution is initially dominated by the first terms in equations (jSJ) and 
(jHJ) so that E' oc M' gas at first, giving a' = 1 before dissipation becomes important. In this 
regime, disk turbulence is fast and driven only by gas accretion. 

Physical quantities can be determined from the normalized variables if we specify the 
corresponding galactic environment. Typical high-redshift star- forming galaxies are charac- 
terized by R w 10 kpc and V w 220 km s" 1 , which give M w 9 x 10 10 M Q and R/V = 45 
Myr. If their accretion rates are 50% of the mass M in one orbit time, 2irR/V = 280 Myr, 
then £ = S/S' ~ 160 M & yr _1 and A = l/(47r) = 0.08. The efficiency of conversion of infall 
energy into disk turbulent energy is not well known. Klessen & Hennebelle (2009) suggest e 
is a few percent to 10%. If we adopt e = 0.1, then t/t' = 90 Myrs, M gas /M^ as = 1.4 x 10 10 M Q , 
E/E' = 1.4 x 10 57 erg, and a /a' = 70 km s -1 . 



3. Results 

The top panels of Figure [1] show solutions to the normalized equations as functions of 
normalized time for star formation efficiency Sq = 0.01. The turbulent energy in the disk 
increases linearly with time at first and then decreases when the dissipation rate becomes 
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larger than the energy accretion rate at about t' — 1. The normalized velocity dispersion in 
the disk starts at a value of unity and decreases when dissipation becomes important. The 
epoch of significant accretion- driven turbulence is t' < 2. 

The unstable mass in a turbulent disk is the Jeans mass, Mj eans ~ cr 4 /G 2 £, which in 
dimensionless units is Mj eans = 1.5tt (a') 4 /M' gas written above. The ratio of the Jeans mass 
to the gas mass is 

M jeans /M gas = (M'j ea jM' g Je 15 /A. (10) 

Initially Mj eans >> M gas and no gravitational instabilities are possible. When Mj eans < 
0.5M gas , instabilities can begin; i.e., each half of the disk can clump into a separate cloud. 
This condition corresponds to M' Jeans /M' gas < 0.5v4/e L5 . For A = 0.08 and e = 0.1 as above, 
this becomes M' 3eans /M gas < 1.26. In Figure [TJ this occurs at t! >0.94 or t > 85 Myr with 
these A and e. 

The time evolution of the normalized instability parameter Q' is shown on the top right 
in Figure [TJ It starts high and decreases as a' decreases and M' gas increases. If significant 
instabilities occur when Q < 1, then this corresponds to Q' < A°' 5 /e " 25 . With A = 0.08 and 
e = 0.1, this requires Q' < 0.5, and it occurs in Figured] at t' > 1.64, or t > 148 Myr. 

The number of giant clouds that form from the instability equals about M gas /Mj eans . 
Note that M gas /Mj eans exceeds the minimum likely value, ~ 2, when the disk is still stable, i.e. 
when Q > 1, for these A and e. Thus as soon as the disk becomes unstable, Q < 1, there is 
enough mass in it to form clouds. The condition for instability is therefore the Toomre Q < 1 
condition, rather than the Jeans minimum- mass condition, M gas /Mj eans > 2. For most A and 
e, this is the case. Figure [2] shows with the top right curve the values of A versus e for which 
M gas /Mj eans = 2 and Q = 1 occur at the same time. A star formation efficiency S = 0.01 is 
assumed. Values of A and e below and to the left of this curve have M gas /M Jeans > 2 when 
the disk is still stable. Thus all of this lower-left region in A — e space begins star formation 
because of the Toomre condition, rather than the Jeans condition. Above and to the right 
of the M gas /Mj eans = 2 curve, the Jeans condition determines the onset of star formation. 
The other black curves in Figure [2] correspond to values of M gas /M Jeans when the first point 
of instability occurs, at Q — 1, also for So = 0.01. Loosely interpreted, these curves give 
the number of giant clouds that form in the disk. The red dashed curves are for a higher 
efficiency, So = 0.1, which does not change the values much. 

The top right panel of Figure [TJ shows the product of the time and the instability growth 
rate. The product u't' increases and exceeds 1 when t' > 0.93. It is generally large by the 
time Q < 1. Note that u't' = ut, independent of A and e, so that when u't' > 1, the 
unstable growth occurs quickly compared to the age of the galaxy. Such quick growth is 
generally applicable in our models. The value of u't' at the onset of instability (Q = 1) 
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is shown in Figure [3j The solid black curve uses So = 0.01 and the dashed red curve uses 
S = 0.1. Both assume Q decreases continuously with turbulent dissipation, i.e., without a 
Q — 1 floor. What is plotted is u't' versus A°' 5 /e ' 25 , which is the value of Q' at Q — 1; 
this is a normalized curve, independent of A and e. For A - 5 /e - 25 < 1.31 in the So = 0.01 
case, u't' > 1 at the onset of instabilities. For fiducial A = 0.08 and e = 0.1, A°- 5 /e - 25 = 0.5 
and u't' = 4.65, as indicated in the figure by the arrow. The rapid turnup in the So = 0.1 
case is because Q' levels off faster at late times than in the So = 0.01 result of the 

decrease in M gas that follows from the higher star formation rate with So = 0.1. 

The top center panel in Figure [1] shows that the normalized star formation rate ap- 
proximately equals the gas accretion rate (= 1 for our normalized quantities) after the time 
t' ~ 3, at which point the ratio of the turbulent dispersion speed to the orbit speed has 
decreased to ~ 0.1 in this So = 0.01 case. Star formation increases so rapidly after t' ~ 2 
that the stellar disk mass soon exceeds the gas mass, which happens at t' = 5.5 in Figure 
[TJ Then the disk enters a near-steady state with a powerful starburst lasting as long as the 
high accretion continues. Disk self-gravity should contribute to the turbulent energy after 
t' ~ 1, when disk instabilities become active, and star formation should contribute to the 
turbulent energy after t' ~ 2, when the starburst begins. 

A problem with these solutions is that a becomes small before the starburst begins 
and before M gas ~ M stars . Observations suggest the opposite, that o is a large fraction of 
the rotation speed in bursting galaxies with M gas ~ M stars . This problem occurs for pure 
accretion-driven turbulence, which always dies out after t' ~ 1. In a real galaxy, gravitational 
instabilities and star formation drive turbulence and prevent Q from dropping much lower 
than 1 (Burkert et al. 2009; Bournaud et al. 2009; Dekel, Sari, & Cervino 2009). We 
simulate this here by allowing Q to decrease as a decreases and M gas increases, but when 
Q = 1, we stop the decrease in a and set a = 7tGE/k, which keeps Q = 1. This means that 
a decreases at first because of the dissipation of accretion energy, but then it increases in 
direct proportion to the gas mass in order to keep Q — 1. 

The bottom panels of Figure [1] show the same variables again as in the top panel, but 
now with this Q > 1 constraint, and also plotted with explicit evaluation of the A and e 
dependencies, so that the variables are physical with normalization to M , V, and R. The 
curve for a/V shows the effect just discussed: the dispersion decreases at first and then 
increases once Q reaches 1 in order to keep Q = 1, as indicated in the right panel where 
Q itself is plotted. The energy decreases at first, from dissipation, but then increases from 
the assumed gravitational instabilities. While this situation is closer to reality than an ever- 
cooling disk, the high a that results keeps the star formation rate low. Then the stellar 
mass builds up very slowly and is always much less than the gas mass over the timespan 
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plotted. Note that SR/(MqV) in the figure barely gets above 0.01, whereas the accretion 
rate in these same units, £R/(M V) = A is 1/(4-71") = 0.08. Thus the star formation rate has 
not yet equilibrated to the accretion rate in this time. Evidently, forcing Q > 1 makes a/V 
more reasonable for high redshift disks, but then star formation is too slow and the stellar 
mass does not readily build up to equal the gas mass. 

To fix this second problem, we have to increase the efficiency of star formation, So- 
Figure H] shows the same normalized variables in the top panel as in Figure [TJ again without 
Q restrictions, and the same physical variables in the bottom panel as in Figure [TJ again 
with Q > 1, but this time So is 10 times larger, 0.1. Now the star formation rate gets large 
quickly, by t' ~ 2, it saturates to the gas accretion rate even in the high-a case (bottom 
panels), and the stellar mass becomes equal to the gas mass within the plotted physical time 
range, at tV/R = 11.56. 

The time when M gas = M star depends on the star formation efficiency. The interesting 
case is when Q has a floor value, so we consider that now. Figure [5] shows the physical 
time, tV/R, when M gas = M star , versus the star formation efficiency, S , on the left and 
the a/V ratio at this time versus So on the right. The solid curves are for a floor value of 
Q — 1 and the dotted curves are for a floor value of Q = 0.7 (beyond which a increases 
proportional to M gas , as discussed above). The different colors are for different turbulence- 
driving efficiencies, e. Both tV/R and a/V decrease with higher So, and they diverge for low 
S . The divergence at low S is because the star formation rate per unit gas mass is so slow 
with Q at its floor value that the accretion keeps the gas mass larger than the stellar mass, 
and our M gas = M star condition for this figure cannot be satisfied. Also at low So, when 
the time for M gas = M star is large, the total gas mass becomes large for the fixed accretion 
rate and a/V also has to be large to keep Q constant. There is a sharp turnaround in a/V 
for So > 0.4 (depending on e) with the Q > 1 curves. In this region, the star formation 
efficiency is so high and M gas = M stars so early, that the turbulence is still from gas accretion 
and Q has not reached its floor value of unity yet. Increasing So beyond 0.4 increases a/V 
because the time gets shorter in the accretion-dominated turbulent regime. These Q > 1 
curves blend smoothly with the Q > 0.7 curves at S > 0.4 (depending on e) because the Q 
floor is not reached in either case for large So and then the Q floor does not matter. 

4. Time-increasing Potential 

The basic equations are modified slightly if the galaxy potential builds up with time 
along with the gas. If we suppose that the potential is V 2 (t/to) instead of V 2 for constant 
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V, then equation [5] becomes 



dE' 



t' - M' disk (E'M' gas ) 



0.5 



S'E'/M' 



(11) 



dt 



gas 



where 



A o.5 V 



(12) 



is the same factor that normalizes t' compared to t (f = t/t ). All of the other equations are 
the same because they do not involve the potential directly, except for Q, which becomes 



because the square root of the potential appears in the epicyclic frequency, k. 

Two solutions of these and the other equations given above are shown in Figure El which 
should be compared with Figure [TJ One solution, with the solid curves, is for a growth of 
the potential until t' — 1 , at which point the galaxy mass is the same as in Figure [TJ and 
then after this the galaxy mass remains constant (i.e., the former equations are used again). 
The second, shown by dashed curves in Figure El is for a continuous growth up to arbitrarily 
large galaxy mass, which uses these new equations for all t'. The first case is a lot like the 
previous solution except for t' < 1, where now the velocity dispersion, a', starts at zero 
instead of ~ 0.8. As a result, the Jeans mass is very small throughout this solution and the 
instability parameter Q' also starts small. The other changes are minor, especially for the 
bottom row of both figures, which assumes Q has a minimum value of 1 and plots unprimed 
quantities with A = 0.08 and e = 0.1. 

In the second case (dashed lines) where the potential continues to grow, the normalized 
velocity dispersion stays large for a long time (see top left panel of Fig. [6]). The rotation 
speed is large too at these late times, higher than before by the factor (t') 0,5 . The increase 
in a' causes the star formation rate to increase more slowly than for a constant potential, 
and the gas mass then gets higher before it equals the stellar mass. The other changes are 
minor for this second case. 

The changes to Figure H] for a time-changing potential (not shown) are qualitatively the 
same as the changes to Figure [TJ that are shown in Figure El The velocity dispersion and Q' 
start low and the Jeans mass is low throughout. The cases when Q has a minimum value of 
1 are hardly changed at all. 

Figure El which showed the time and relative velocity dispersion when M gas = M star 
in the Q > 1 case, is also hardly changed when the potential increases with time. This 
is characteristic for Q > 1. However, the ratio of a to the full rotation speed is now 



Q' = (3t')°V/M' 



gas 



(13) 
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a I (Vft/to] 1 ^ 2 ) instead of a/V. This change can be significant, as seen in Figure [71 The 
left-hand panel is for an increasing potential up to t — t and then a constant V. The right 
hand panel is for a continuously increasing potential, with no limit. The curves in the left 
panel are almost the same as in Figure because the time when Sq < 0.3 is larger than 
to and then the galaxy mass is the same as before. Recall that to — e°' 25 i?/ (A ' 5 V), so 
t/t = (A°- 5 /e°- 25 )(tV/R) ~ 0.50(tV/R) for A = l/Air and e = 0.1. This means that the time 
tV/R plotted in the left hand panel of Figure [5] should be multiplied by 0.5 to convert to 
t/t . The curves in the right panel have lower values of a/V because the mass is larger at 
larger times, and so the rotation speed is larger, when So is small. The implication of this 
change is that observed ratios of dispersion to rotation speed for galaxies with M gas = M stars 
permit smaller star formation efficiencies if the galaxy mass continuously increases during 
gas accretion. The primary reason is that relative dispersions of several tenths can occur 
later in galaxies with M gas = M stars if the galaxy mass is larger. 



5. Comparison with Observations 

High redshift galaxies are not yet observed at the phase that is most relevant to this 
paper, namely, when turbulence is still driven by gas accretion and star formation has not 
yet begun. Selection effects limit the observation of galaxies to the starburst phase, at which 
point a significant amount of turbulence should be driven by the dynamics that triggers the 
star formation, most likely gravitational instabilities, and also by the young stars themselves. 
Thus the fast turbulence observed in young galaxies so far is probably not from accretion 
energy, but from gravitational instabilities and young stars. However, the turbulence that 
was present before the giant star-forming clumps form, which caused the disk Jeans mass 
to be so large and gave these clumps their enormous masses, could have been accretion 
energy. We predict that observations of extremely young, gas-dominated disks will show 
high turbulent speeds even before star formation begins. 

In the simple model presented here, accretion energy causes the turbulent motions 
during the first few rotation times for a steady, high accretion rate. In general terms, 
this phase lasts for t acc = 2e°' 25 R/(A°' 5 V) considering that t' ~ 2 at this time (from Figs. 
CDandlU). Writing A = £R/(VM Q ), the duration becomes t acc = 2 {e°- 5 RM /[£V}) - 5 = 
2 (e°- r °R 2 V/^G}) - 5 = 2.2 (e a5 GM 2 /[^ 3 ]) a5 for accretion rate £ in physical units. Thus the 
accretion-driven turbulence phase lasts longer for higher efficiency e and lower accretion rate, 
and for more massive or larger galaxies at a given rotation speed. For the typical parameter 
values discussed elsewhere, t 

acc rsJ 

180 Myr. 

The normalized time when the disk becomes gravitationally unstable, i.e., Q — 1, which 
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means Q' = A°- 5 /e - 25 = 0.5 for A = 0.08 and e = 0.1, is found to be t' ~ 1.7 (1.64 in Fig. □ 
with Sq = 0.01 and 1.74 in Fig. HJwith Sq = 0.1). These are for the cases with no floor in Q. 
The normalized gas mass at this time is M' gas ~ 1.6, (1.62 for 5*0 = 0.01 and 1.53 for So = 0.1) 
which means the physical gas mass is M gas = M gas A°- 5 e°- 25 M = 0.25M = 2.3 x 10 10 M 
for this normalization. This is the gas mass at the beginning of the starburst phase. By 
the time the gas and stellar masses are equal, which is appropriate for the current gas 
observations (e.g., Tacconi et al. 2009), t' — 3 to 5 (5.5 for S = 0.01 and 3.12 for S = 5.5) 
and M' gas = 2.7 to 1.6, respectively, making M gas ~ 3.9 to 2.3 x 10 10 M Q . With a Q = 1 
floor and Sq = 0.1 (bottom panel of Fig. HJ), tV/R = 11.6 when M gas = M sta rs, and then 
M gas = 0.46M = 4.1 x 10 10 M . 

These gas masses are appropriate for the start of the star formation phase, between the 
time of first instabilities and the time when the stellar mass has built up to be comparable to 
the gaseous mass. The masses are large in our model because the young disk is stabilized by 
turbulence that is driven by accretion energy. This energy source explains how the gas mass 
can build up to such large values without first turning into stars. The youngest disks require 
stability like this or else the gas-to-star ratio will always be low. We obtain the observed high 
ratio and the observed gas masses for reasonable values of the accretion rate and efficiency 
of turbulence driving. We also obtain the high Jeans masses of the first star-forming events, 
which is also an indication of fast turbulence. 

Our conclusion that accretion dominates turbulence only during the very earliest phase 
of galaxy growth may be derived most simply from the dimensionless equations in Sec- 
tion [2j From equation [5], we see that energy gains exceed energy dissipation when 1 > 
M^ isk (i?'M gas ) ' 5 . From equation [3, M^ isk = t'. Also at early times, there are not many 
stars and equation [8] gives M gas ~ t'. Similarly, from the first term in equation [5], E' ~ t'. 
From all of these we derive t' < 1 when turbulence is driven by accretion. This translates 
to a physical time t < t = e°- 25 R/ (A°- 5 V), or, after substituting A = ^R/{VM ), we get 
t 2 < e°- 5 RM / (£V), as derived above. 

Dekel, Sari & Ceverino (2009) consider disk turbulence driven by accretion, disk in- 
stabilities, disk clump stirring, and star formation. For the accretion-dominant phase, they 
derive a timescale in their equation 40, which is ~ M^s^a 2 / (£V 2 ). This is essentially the 
same as we get because we both assume equality between the turbulent energy and the ac- 
cretion energy in this phase. The simple derivation in the previous paragraph eliminates a 
in favor of global galactic quantities because a is what we wish to determine. Most of the 
discussion in Dekel et al. is about disk self-regulation, which is more simply treated in the 
present paper by assuming the instability parameter Q has a minimum value of order unity. 
Self-regulation occurs at a later phase, if > 1, when accretion is not strong enough to drive 
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turbulence in comparison to internal disk processes and star formation. Observations so far 
see primarily this late phase. Lehnert et al (2009) comment on accretion- driven turbulence 
after observing spectra of redshift z ~ 2 galaxies. They note that the observed turbulence is 
too dissipative to be driven by a reasonable accretion rate and conclude that star formation 
dominates instead. Our model agrees with their assessment, as the starburst phase occurs 
much later than the accretion-driven turbulence phase, at t' ~> 3 (Fig. [p. 

Klessen & Hennebelle (2009) suggest that turbulent accretion can power turbulence in 
a wide variety of conditions, including the outer parts of modern galaxies where the star 
formation rate is low. For high redshift galaxies, they conclude that the observed turbulence 
inside individual clumps can be driven by clump gas accretion if the total rate is 10 — 50 M Q 
yr _1 . They consider this reasonable as it is comparable to the observed star formation rate. 
However, they also caution that their model might not apply to the whole disk of a clumpy 
high redshift galaxy, as clump coalescence and minor mergers might dominate accretion and 
star formation in that phase. 



6. Conclusions 

Gas accretion to a young galaxy disk can drive turbulence for a time t acc ~ 2 (e°' 5 GM 2 /£,V 3 ) 
where e is the fraction of the accretion energy going into disk turbulence and £ is the ac- 
cretion rate. After this time, disk turbulence should be driven by gravitational instabilities 
and star formation. The first instabilities should produce only a few giant gas clumps that 
dominate early star formation until the relative gas fraction in the disk decreases. The 
star formation rate equilibrates to the accretion rate in 1 or 2 t acc , depending on the star 
formation efficiency (So). 

Observations of a relative high turbulent speed compared to the rotation speed, and of 
a disk gas mass comparable to the disk stellar mass, require that the galaxy-averaged star 
formation efficiency has to be large compared to the modern galaxy-averaged star formation 
efficiency. This efficiency is defined here as the star formation rate per unit dynamical time 
in the disk, and per unit total galaxy area. Thus a high average value means either that 
the local efficiency is high and the g equals the galaxy area, or the local efficiency is 

more normal and the g cLS 3X6 cl IS smaller than the galaxy area. Our one-zone model cannot 
distinguish between these possibilities. The high efficiency causes the stellar mass to build 
up to the gas mass before the velocity dispersion in a Q = 1 gas disk has dropped below a 
few tenths of the rotation speed. This criterion is easier to satisfy if the total galaxy mass 
increases with the gas mass (Fig. [7]). 
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The high rate of star formation in young galaxy disks is the result of a rapid instability 
compared to the galaxy age and a high gas mass compared to the galaxy mass. The rapid 
instability follows from our simple model (Fig. 3). The instability growth rate, u, compared 
to the rotation rate, V/R, is uR/V ~ (M gas /M) (V/cr). This is large because the relative 
gas mass, (M gas /M), is large for young disks, even though the relative velocity dispersion, 
(cr/V) is somewhat large too. As (cr/V) drops because of turbulent dissipation, and the 
accretion continues, uV/R becomes even larger. The star formation rate declines only when 
the accretion rate declines. 

We conclude that gas accretion is a good source of energy for ISM turbulence in the 
earliest phases of galaxy growth. It is quickly replaced by other sources after a few rotation 
times, and is not likely to be the source of turbulence that is observed at high redshift in 
starbursting systems. Accretion-driven turbulence in a young disk is important because it 
provides the initial stability that allows the gas to build up to a large mass before star for- 
mation begins. It also gives the observed large masses for star-forming clumps. Observations 
of gas-rich disks before the starburst phase should show the predicted high turbulent speeds 
that come from accretion. 

We are grateful to the referee for helpful suggestions. 
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tV/R tV/R tV/R 

Fig. 1. — Normalized solutions to the equations of accretion, dissipation, and star formation 
are on the top, and physical solutions with Q regulated to be larger than or equal to 1 are 
on the bottom, both for the case where the star formation efficiency, Sq, equals 0.01. The 
left panels show energy and velocity dispersion, the middle panels show masses and the star 
formation rate, the right panels show the stability parameter and the product of the growth 
rate and time. 
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Fig. 2. — The curves show the ratio of the gas mass to the turbulent Jeans mass as a function 
of the normalized accretion rate, A = £R/(VM Q ) for physical rate £, and the efficiency e 
for the conversion of accretion energy into turbulence. The solid black curves assume a star 
formation efficiency of Sq = 0.01 and the three red dashed curves assume So = 0.1 for the 
M gas /Mj eans = 2, 5, and 20 cases. This ratio of masses should be about equal to the number 
of giant clumps formed in the first epoch of instabilities. 
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Fig. 3. — The product of the normalized growth rate and the normalized time, which is the 
same as the product of the physical growth rate and time, versus a combination of A (ac- 
cretion rate) and e (turbulence driving efficiency) that appears in the normalized instability 
parameter Q' . For A = 0.08 and e = 0.1, as discussed in the text, this combination equals 
0.5, which is shown by the arrow. The large values of uo't' suggest that the instability should 
be rapid once the disk goes unstable. The solid black curve is for So = 0.01 and the dashed 
red curve is for So = 0.1. 



- 19 - 



> 

"6 




1 


1 




/ / 




// 




// 




CO' t'A 




// 




// 






— 1 1 1 



2 4 6 8 10 12 14 
tV/R 



2 4 6 8 10 12 14 

tV/R 



2 4 6 8 10 12 14 
tV/R 



Fig. 4. — Same as Figure 1, but with higher star formation efficiency, Sq = 0.1. With 
higher efficiency, the stellar mass increases sooner and faster, reaching equilibrium with the 
accretion rate earlier. High So is required to get M gas ~ M sta rs and large a/V, as observed, 
when the velocity dispersion is regulated by Q > 1. 
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Fig. 5. — The relative time, tV/R, (left) and relative velocity dispersion, cr/V, are shown 
versus the star formation efficiency, Sq, at the time when M gas = M stars . This epoch of 
mass equality is considered because it is representative of disk galaxies at z ~ 2. In all of 
these solutions, the Q parameter has a floor value (1 or 0.7), which is enforced by making a 
increase in proportion to M gas . For low efficiencies, the star formation rate cannot keep up 
with the accretion rate, so M stars is never equal to M gas . For high efficiencies, stars build up 
quickly, decreasing the time when a significant mass of stars appears, and decreasing the gas 
mass at that time, from which the lower cr/V follows. The turn-around in cr/V at high So is 
from the occurrence of M gas = M stars before the Q floor is reached; elsewhere, M gas = M stars 
after the Q floor is reached. 
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Normalized Solutions 
SF Efficiency = 0.01 




Physical Solutions, Q > 1 
A = 1/4ti ; e = 0.1 
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Fig. 6. — Models analogous to those in Figure 1 but now with a potential that increases 
linearly with time. The solid curves are for a case where the potential grows only until a 
dimensionless time if = 1, at which point the galaxy mass remains the same as in Figure 1. 
The dashed curves are for a case where the potential continues to grow for all time. The 
primary effect of a growing potential is to decrease the velocity dispersion, Jeans mass, and 
stability parameter Q at small times. The lower panels, which assume Q is regulated to 
remain larger than 1, are not significantly affected by the change. 
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Fig. 7. — The ratio of the dispersion to the instantaneous rotation speed is shown at the 
time when M gas = M stars , as in Figure 5, but now with a time- increasing potential. The solid 
and dashed lines are the same as before (Q > 1 and Q > 0.7 self-regulation, respectively), 
and the three cases for the efficiency of conversion of accretion energy into turbulent energy, 
e, are also the same. On the left, the potential increases only until t' — 1 after which the 
mass remains the same as in Figure 5. On the right, accretion continues for all time. The 
left-hand panel is not changed much from before, but the right-hand panel suggests that at 
low Sq, which corresponds to long times, a/V(t) is lower when the potential continues to 
increases than when the potential was constant (Fig. 5) or stops increasing at t' = 1 (left 
panel). This is because the rotation speed is large when the potential is large. 



